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abstract Results have been obtained on the quasi-elastic spectra of neutrons scattered from pure water, a 20% 
agarose gel (hydration four grams H 2 0 per gram of dry solid) and cysts of the brine shrimp Artemia for hydrations 
between 0.10 and 1.2 grams H z O per gram of dry solids. The spectra were interpreted using a two-component model 
that included contributions from the covalently bonded protons and the hydration water, and a mobile water fraction. 
The mobile fraction was described by a jump-diffusion correlation function for the translation motion and a simple 
diffusive orientational correlation function. The results for the line widths T(Q 2 ) for pure water were in good agreement 
with previous measurements. The agarose results were consistent with NMR measurements that show a slightly 
reduced translational diffusion for the mobile water fraction. The Artemia results show that the translational diffusion 
coefficient of the mobile water fraction was greatly reduced from that of pure water. The line width was determined 
mainly by the rotational motion, which was also substantially reduced from the pure water value as determined from 
dielectric relaxation studies. The translational and rotational diffusion parameters were consistent with the NMR 
measurements of diffusion and relaxation. Values for the hydration fraction and the mean square thermal displacement 
(u 2 ) as determined from the (7-dependence of the line areas were also obtained. 

INTRODUCTION 

Water is the most abundant molecule in biological systems, 
ranging from -50% to 90% by weight in most organisms, 
tissues, and cells, depending on the physiological state. The 
role of water in biological systems has been a subject of 
controversy; however, its role in biochemical and biological 
processes is certainly one of the most important unsolved 
problems in biology (1). The controversy arises in part 
from the lack of fundamental knowledge of the micro¬ 
scopic structure of pure water and water in biological 
systems. It also arises from the fact that the results of 
measurements with the various technologies used to study 
the properties of water in biological systems have not been 
consistently integrated with one another and with any 
reasonable model of water structure and dynamics. 
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The physical models used to describe the behavior of 
water in biological systems range from a simple bound-free 
exchange model (2-4) to the polarized multilayer hypothe¬ 
sis of Ling (5-10). Although experiments have been per¬ 
formed to test these views, no general consensus has 
emerged. In this paper, we review the NMR experiments 
on biological water, compare their results with those 
obtained from neutron scattering, and show that neutron 
scattering gives new results not obtainable by other meth¬ 
ods. 

MEASURING METHODS 
Neutron Scattering Technique 

The technique of quasi-elastic neutron scattering (QNS) is a powerful 
method for studying the nonperiodic (diffusive) motion of atoms in solids 
and liquids (11,12). The results of QNS spectra measurements are 
usually interpreted within the framework of the Van Hove theory 
(13, 14), which relates the scattering law S(Q y a>) to the space-time 
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correlation function G(r,t) of the scattering particles. This correlation 
function can be calculated for various models for the assumed motion of 
the diffusing particles (e.g.. Brownian, jump-diffusion, and oscillatory 
diffusion). The shape of the expected QNS spectra can be determined 
from the correlation function, and the parameters of the models deter¬ 
mined by a fit to the experimental spectra. These parameters can then be 
related to the microscopic properties of the environment of the scattering 
atoms, (e.g., residence time, jump length, and rotational diffusion coeffi¬ 
cient). 

This method holds special promise for the study of the properties of 
water in biological systems. The structure of pure water can be studied 
with x-rays (15,16), but these methods cannot be easily extended to 
noncrystalline biological systems, due to their heterogeneous molecular 
composition. Nuclear magnetic resonance (NMR) techniques have been 
extensively used to study the properties of water by measurement of the 
various relaxation times and the translational diffusion coefficient (3). 
The NMR relaxation parameters are determined by the environment of 
the proton moments (dipole-dipole interactions with other protons, 
located on other water molecules or within the cellular macromolecules) 
(17, 18), and the diffusive properties are determined by the intermolecu- 
lar interactions as well as by the interaction with the biomolecular 
structures which can form obstructions and barriers for diffusion (19). 

NMR measurements on the properties of water in biological systems 
have not yet provided a clear picture of the microscopic disposition of 
water within the cell. This is due to the heterogeneous character of the cell 
as well as to the fact that a single NMR measurement is carried out over a 
period of time, t, of the order of a few milliseconds or more (20-22). 
During this measuring time, the water molecules are diffusing within the 
cell, and their NMR parameters are determined by averages over 
distances y/2Dt =* 1 fim where D is the diffusion coefficient of water. An 
interpretation of the results requires a detailed modeling of the heteroge¬ 
neous environment of the cell, and there has not been universal agreement 
on the essentail features of such a model. 

The QNS method gives information on the properties of the scattering 
atoms that depends directly on the microscopic environment of these 
atoms. The correlation function for the atomic motion is explored over a 
space domain on the order of a few Angstroms and for f~10~ 12 s. The 
positional correlation of the atom determines the scattering. This is to be 
contrasted with the NMR measurements that give information on the 
correlation function for the orientation of the protonic moments (23). 

In this paper, we present QNS spectra for water, an agarose gel, and 
four hydrations of Artemia cysts. The spectra are composed of two 
components, one that is elastic within the spectrometer resolution and 
another that is broadened by the diffusive motion. These spectra have 
been interpreted with diffusive correlation functions for the translational 
and rotational motion of the mobile water, and hydration (or “bound”) 
fractions have been determined from an analysis of the line areas. 

Sample Preparation 

The water used for the preparation of all samples was double distilled and 
deionized. The pure water sample was contained in an aluminum sample 
chamber. 

Agarose gel is a polymer (C 12 H 18 0 9 ) n of known structure (24). It 
contains 4 OH groups and 14 covalently bonded protons in each repeat 
unit. It was selected for a pilot study to determine if the spectrometer 
resolution was sufficient to separate the quasi-elastic line due to water 
protons from the elastic line of the polymer protons. In addition, NMR 
data on the diffusion coefficient and relaxation times have been published 
(25). Gels of 20% agarose by weight in water and deuterium oxide were 
prepared as follows. The appropriate amounts of de-ionized, distilled 
water or deuterium oxide, and electrofluorescence grade agarose powder 
were weighed, and then mixed in a small beaker. The sample was then 
covered and heated to 90°C in a water bath and maintained between 85° 
and 95°C for 30-60 min. The samples were then poured into preheated 
aluminum sample holders and allowed to set at room temperature for a 
minimum of 12 h before placing them in the spectrometer. 


The cysts of Artemia , a primitive crustacean called the brine shrimp, 
consist of an inner mass of ~4,000 eucaryotic cells surrounded by a 
complex noncellular shell. Their ultrastructure, biochemistry, and physi¬ 
cal properties have been described in some detail (26), and they have been 
used for the study of a wide variety of biological problems. Their utility 
for studies on cellular water arises from the ability of these cysts to 
undergo virtually complete dehydration in a reversible way (27) and to 
tolerate the conditions necessary to carry out the QNS study (28). The 
cysts used in these experiments were purchased from San Francisco Bay 
Brand, Inc. (Menlo Park, CA.) Approximately 90 ± 3% of this popula¬ 
tion, after processing, produced viable larvae when incubated at 25°C in 
seawater for 72 h. Large quantities of cysts were prepared by a thorough 
washing procedure carried out at 2-4°C, to suppress metabolic activity 
(29). The cysts were then dried at room temperature, separated by size 
(average diameter ~200 /*m) and stored in a desiccator over CaS0 4 until 
used. 

Various water contents (hydrations) were achieved by allowing the 
cysts to equilibrate with NaCl solutions of different concentration. The 
cysts were immersed overnight at 0°C. Because the inner shell of the cyst 
is impermeable to NaCl (28), the amount of water taken up is determined 
strictly by water activity of the NaCl solution. By this procedure, 
hydrations between 0.1 and 1.2 grams of H 2 0 per gram of dry cyst (g/g) 
were obtained. Maximum cyst hydrations are ~1.4 g/g, at which the 
cellular water content is ~ 1.7 g/g (~63% by weight). Details of hydration 
methods have been published (29). 

The cysts were removed from solution and the surface water removed 
quickly. They were allowed to equilibrate with each other for 30-60 min 
in a closed vial; then they were transferred to an aluminum sample 
chamber for QNS measurements. A sample was also taken for prelimi¬ 
nary determination of their water content by gravimetric procedures. 
(Drying was carried out at 100°-105°C for at least 6 h.) A final value of 
the water content was obtained from measurements on cysts taken from 
the sample chamber after the QNS experiment was completed. The two 
hydration values were always in close agreement, showing that cyst 
hydration was constant during the course of the experiment. Hatching 
tests showed that the viability before and after neutron exposure was the 
same as for controls that had not been exposed. Thus, exposure to 
neutrons had no effect on viability, and the results can be regarded as 
those for living intact cells. 

Neutron Spectrometer 

All of the spectra reported in this paper were taken on a triple-axis 
spectrometer (labeled HB2) at the High Flux Isotope Reactor (HFIR) at 
Oak Ridge National Laboratory. In these experiments, the monochroma¬ 
tor and analyzer were pyrolytic graphite [(002) reflection]. Neutrons 
from higher-order reflections were removed with a beryllium filter. The 
analyzer was oriented to accept neutrons with an energy of 0.95 THz 
(~3.9 meV), and the neutron energy transfer was scanned by scanning the 
diffraction angle of the monochromator. Collimators with an angular full 
width at half-maximum of 20 min of arc were placed before and after the 
sample. The resolution function was determined from a run on a 
vanadium sample and was Gaussian with a full width at half maximum of 
0.025 THz. The spectrometer scans were controlled by a PDP-8 com¬ 
puter. (Digital Equipment Corp., Marlboro, MA). 

The intensity of the scattered beam was determined with a 3 He 
detector, interfaced with the computer. The counting time was deter¬ 
mined by a low efficiency (uranium) fission detector that monitored the 
intensity of the incident beam. The data consisted of a print-out of the 
number of scattered neutrons obtained in the 3 He detector at each value of 
momentum and energy transfer during the time required to obtain a 
preselected number of monitor counts. 

The sample holders were made from a standard aluminum alloy 
(2024), and they were machined so as to provide a sample chamber either 
1 or 2 mm in thickness, and with lateral dimensions ~4 cm x 4 cm. The 
front and back covers had a thickness of 1 mm, and the chambers were 
sealed with indium wire gaskets and aluminum nuts and bolts. The beam 
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(-2.5 x 2.5 cm) always passed symmetrically through the sample 
chamber. 

The thickness of the sample chambers was 1 mm for the pure water 
sample and for the agarose-H 2 0 gel. It was 2 mm for the agarose-D 2 0 gel 
and the brine shrimp cysts. These thicknesses were chosen so as to make 
the thickness, measured as a fraction of a neutron mean free path, roughly 
equal to 1 /2 for all samples, which gives a transmission of -60%. (The 
mean free path of neutrons in pure water is -2 mm.) 

The wavelength of the 0.95 THz neutrons used in these experiments 
was not sufficiently long to avoid all coherent reflections from aluminum. 
Reflection from the (ill) planes was possible at a Bragg angle of 6 = 
79°. The fraction of the beam scattered from this reflection was estimated 
to be -3.5%, which would then have to be scattered by the water at rather 
large 0-values to contribute to the QNS spectra. We estimate that these 
neutrons make a small contribution to the background, but do not 
otherwise affect the spectra. All data were collected with the samples at 
the temperature of the spectrometer bay, which was T = 21° ± 1°C. 

DATA 

Because the thicknesses of the various samples were 
adjusted to be roughly one-half a neutron mean free path, 
the number of counts from the different samples for each 
momentum transfer hQ and energy transfer ho) were 
roughly the same for a given monitor setting. Approxi¬ 
mately 100-400 counts per point were obtained in the far 
wings of the spectra, where the background and inelastic 
scattering are dominant, and -1,000-3,000 counts were 
collected at the peak of the quasi-elastic line. 

The spectral scans were usually taken at five Q values: 
Q = 0.7, 1.0, 1.2, 1.5, and 1.9 A -1 and for —40 values of 
energy transfer hco between ±0.200 THz. One scan was 
taken of the Artemia cysts of hydration 1.2 g/g at Q = 0.5 
A -1 . The separation of points near the quasi-elastic peaks 
was Aco = 0.005 THz. The samples from which spectra 
were obtained and the hydration of these samples are 
presented in Table I. Spectra were also obtained on an 
agarose-D 2 0 gel whose hydration level was adjusted so 
that the agarose content in mole percent was equal to that 
of the agarose-H 2 0 gel. 

Typical raw data for the H 2 0 and D z O agarose gels as 
obtained from the computer print-out are shown in Fig. 1. 
The number of counts as a function of energy transfer for 
0- 1.2A- 1 is shown. The broadening for the H 2 0 gel that 
is due to the mobile water is evident. 


TABLE i 

WATER CONTENT OF NEUTRON 
SCATTERING SAMPLES 


Sample 

Hydration 

gm H 2 0/gm dry solids 

h 2 o 

CO 

Agarose gel 

4.0 

Brine shrimp cysts 

1.2 


0.78 


0.31 


0.10 


ANALYSIS AND RESULTS 

Models for Pure Water 

In liquids such as water where atoms are associated with 
nearest neighbors by hydrogen bonds, the translational 
motion of the atoms is often described by a “jump- 
diffusion” model, which gives a Lorentzian scattering law 
(11,12): 


St«2,«) 


r t (fi) exp(-e 2 <^» 




where hQ is the momentum change of the scattered 
neutron, T t (0) is the full-width at half-maximum, and 
exp( — ^(Wm 2 )) Is the Debye-Waller factor. One form of 
this model assumes that the molecules perform oscillations 
at quasi-equilibrium positions for a time r, and that they 
then diffuse to a new position during the time r,. If r, « r, 
then 


r t (0) = 2{0 2 Z)+ [1 -exp(-e 2 <^»]/r}/(l +Q 2 Dt). (2) 

The mean square jump length is (f 2 ) == 6 Z>,t„ £>, is the 
diffusion coefficient during r h and D = (1/6)((/ 2 )/t) + 
( u m ) / T is the macroscopic diffusion coefficient. 

A modification of this model (30, 31) in which the 
diffusion during time r, is replaced by an instantaneous 
jump to a new site at a distance /, and the oscillations 
treated with the Gaussian approximation, gives the width 

T t = 2Q 2 D/(\ + 0 2 Z>t), (3) 

provided that the probability /?(/) for a jump of length / is 
p(/) a / exp( —/// 0 ). The macroscopic diffusion coefficient 
is D = I q 2 /t = </ 2 >/6t. 

Eqs. 2 and 3 can both be fitted to the data over the Q 
range that can be easily studied, but there are some reasons 
to prefer Eq. 3 (30, 31). Because Eq. 3 has fewer parame¬ 
ters, we have adopted Eq. 3 for the interpretation of our 
data. 
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Figure 1 Neutron spectra at Q — 1.2 A 1 for a 4 g/g gel of agarose, o, 
D 2 0 gel; A, H 2 0 gel. The spectrometer resolution function is Gaussian, 
with width shown. 
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The jump-diffusion model does not take the rotational 
motion of the water molecule into account, and it may 
seem surprising that its use leads to diffusive parameters in 
good agreement with measurements of the translational 
diffusion by other methods (e.g., tracer, NMR). Springer 
suggests (reference 11, p. 46) that the rotational motion in 
bulk water has a negligible effect on the line widths, since 
the reorientational correlation time is several times larger 
than the residence time r of the jump-diffusion model, so 
the line width is determined principally by the translational 
motion. (This will be considered further in the Discussion 
section.) Irish et al. (32) have shown that the rotational 
motion makes a significant contribution to the line width at 
1°C but does not contribute to the width at 22°C. 

To assess the importance of the rotational motion to our 
measurements, we have adopted a scattering law for the 
rotational contribution based on the following assumptions: 
(a) The rotational and translational motions are uncorre¬ 
lated, so that the Gaussian approximation holds (reference 
14, pp. 375 ff.) ( b) The scattering law £ t ((J,a?) for the 
translational motion is given by the jump-diffusion model, 
Eqs. 1 and 3. (c) The scattering law S r (Q 9 w) for the 
rotational motion is given by a rotational diffusion model in 
which the orientation of the molecular axes undergoes 
Brownian motion described by a rotational diffusion coeffi¬ 
cient D r (reference 11, pp. 66, 67): 

S r (Q,<*)-jHQfl)n<*) 

+ 3 j 2 ,(Qa) — Ff/ ^ + (higher terms) (4) 

"• + (t) 


corresponds to water of hydration, which is tightly bound 
to macromolecules and structures within these systems, 
then we expect the protons in this phase to scatter nearly 
elastically. For convenience, we will henceforth refer to 
water whose quasi-elastic line width is small compared 
with the spectrometer resolution as bound water. This 
means only that its diffusional motion is ^4 times slower 
than pure water. This bound fraction may, however, be 
undergoing some high-frequency oscillations as in a solid, 
so there will be a Debye-Waller Q dependence to the 
intensity of this component. We assume that this component 
has the scattering law: S b ((?,<o) = exp(-0 2 ( w b 2 ))5(<o). 

The second component is assumed to be mobile and to 
diffuse with a scattering law given by Eq. 5. Because this 
phase may still be affected in some way by the presence of 
the other components of the system, we will not assume 
that the free parameters in Eq. 5 are the same as those 
found for pure water. If p is the fraction of the water that is 
in the bound phase, then the total scattering law is given by 
S(Qv) =pS b (Q, o>) + (1 - p)S n {Q, o>). 

The width of the quasi-elastic line given by Eq. 5 is on 
the order of the spectrometer resolution. Therefore, the 
observed partial differential cross section will be the convo¬ 
lution of the “ideal” partial differential cross section with 
the resolution function for the spectrometer. Because the 
measured resolution function is Gaussian, the observed line 
shape will be a Gaussian plus two Gaussian-Lorentzian 
convolutions: 

<*>) = A h G(u>) + A m G(a>)* [T t (o?) + T r+t (o)] 

where 


where j n (z) is a spherical Bessel function of the first kind; a 
is the molecular “radius”; T r = 4 D r ; 6(o;) is the Dirac delta 
function; and the higher terms will be neglected. 

The resulting scattering law S m (Q,w) is the convolution 
of S T with S t : 


S m (Qy <*>) - 


1 j 2 o(Qa)T t c\p(-Q 2 (u 2 m )) 3 , 2/r% 

2* a ; 2 + (r,/2 ) 2 + 2 ^ J1 (Qa) 

(r r + r,)exp(-g 2 <H 2 )) 


*■♦ 


( 5 ) 


Bound-Free Model 

Because the agarose and Artemia cyst systems are hetero¬ 
geneous, the QNS spectra from them cannot be expected 
to fit a single component model. Instead, we expect the 
protons from the various cellular environments (macromo- 
lecular, hydration, and mobile water) to have their own 
scattering laws that contribute to the total observed partial 
differential cross section. 

The simplest model that takes into account the hetero¬ 
geneity of these complex systems is one in which there are 
two phases of water. If we assume that one of these phases 


G( a?) = exp [—(a>/a) 2 ], 

Z, t (a>) = 1/(1 +(*/r t ) 2 
£,.,(«) = i/d + [»/(r r + r,)] 2 | 


and 

G(o?)*£(a>) = J G(o/)L(a> — a?')da?'. 

Results for Pure Water 

The number of counts as a function of w at each value of Q 
were fitted with a function given by F(w) = aw + b + A • 
G(o>) * L t (a?) where aw 4- b represents the background. 
(The rotational contribution to the line was neglected; this 
will be taken up in the Discussion section.) The fitting was 
performed with a generalized least-squares search routine. 
The algorithm is Marquardt’s search technique (33), 
which combines the best aspects of the Taylor’s series 
(Newton-Raphson) method and the gradient search 
method (34). Four parameters were varied: a, b, A , and the 
width T of L x (w). Values of x 2 per degree of freedom were 
in the range 0.85 to 1.05. The area of the Lorentzian 
component followed a Debye-Waller dependence as 
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described in Eq. 1. The mean-square thermal displacement 
(u m 2 ) was 0.33 A 2 . These data are plotted in Fig. 2. This 
value of (u m 2 ) is in good agreement with the results of 
other workers (35). (Some investigators [see reference 12, 
p. 349] prefer to call 6 (u m 2 ) the mean-square thermal 
displacement.) It should be noted that the range of values 
that have been reported for ( u m 2 ) is quite large (reference 
12, p. 349), and there is not always good agreement 
between QNS results and other methods. 

The full width at half maximum T is plotted in Fig. 3. 
The data were fit with the model given by Eqs. 1 and 3, 
with the results D = 2.4 x 10" 5 cm 2 /s and r = 1.2 x 10~ 12 
s. The fitted line in Fig. 3 is a plot of Eq. 3 with these 
parameters for pure water. The self-diffusion coefficient 
(/)) is in good agreement with the results of other workers 
using QNS and other methods (11, 12, 36). The value of r 

10~ 12 s agrees with other QNS results (11, 12). 

Results for Agarose 

The agarose polymer contains a significant number of 
chemically bonded protons. To separate the contribution of 
the agarose protons from the water, one can take advan¬ 
tage of the appreciably lower cross section of deuterium 
oxide compared with water. As described in the Measuring 
Method section, gels of agarose in water and agarose in 
deuterium oxide were prepared. The latter gel was used to 
determine the contribution of the agarose protons to the 
QNS spectrum. 

The deuterium oxide gel data could also be fitted with a 
Gaussian plus a Gaussian-Lorentzian convolution, where 
the Gaussian had the width of the spectrometer resolution. 
The full fitting function was of the form F(co) = av + b + 
A • G(cjj) + B • G(o)) * L b (co). The amplitudes (A and 
B) and the width of the Lorentzian L B fully described the 
scattering from the agarose protons in the gels. The area of 
the Lorentzian contribution was of the same order as that 
of the Gaussian, and its presence was probably due to the 
motion of the OH side groups of the molecule. 

The data obtained from the D 2 0 gel were used to 



Figure 2 Debye-Waller plot for pure H 2 0: area vs. Q 1 . The slope gives 
the mean-square radius of the “thermal cloud” for the oscillatory motion 
of the protons. 



Figure 3 Line width T (in millielectron volts) vs. g 2 (in reciprocal 
Angstroms squared) for pure water, the 4 g/g agarose-H 2 0 gel and the 
1.2 g/g/ Artemia cysts. Water data for g 2 < 2 A~ 2 have been omitted. 
The curves drawn through the data are obtained from diffusion models 
and parameters described in the text. 

remove the contribution of the agarose protons from the 
H 2 0 gel data. The parameters (amplitude and width) of 
the Lorentzian component of the D 2 0 gel, properly scaled 
according to agarose content, were used as fixed input 
parameters for the fits to the H 2 0 gel data. The Gaussian 
amplitude of the H 2 0 gel will contain a contribution of 
unknown amount from the bound water as well as a 
contribution from the agarose protons that is known from 
the measurements on the D 2 0 gel. The form of the fitting 
function for the water gels was thus: F(co) = au + b + A • 
G( co) + B • G(u>) * L b ( co) + C • (7(a?) * Lt(oj), where B 
and L b were fixed at values determined by the D 2 0 spectra, 
and the value of A contained contributions from the bound 
water and the agarose protons. The amplitude C of the 
Lorentzian due to the mobile water was -1-5 times B , 
depending on the value of 0. The areas of the various 
contributions are given in Table II. 

To estimate the amount of bound and free water in the 
agarose gel the areas of the Gaussian and Lorentzian 
components at 0 = 0 are used. The area as a function of 0 
will be called 5(0), which is the integral of the scattering 
law 5(0,0?) over o?. At 0 = 0 this area is proportional to the 
number of nuclei contributing to that portion of the 
spectrum (i.e., Gaussian or Lorentzian). If the scattering 
from the bound water had a Lorentzian contribution, as 
does the scattering from the agarose protons, these areas 
would not be strictly proportional to the corresponding 
number of nuclei. We do not have a good estimate of the 
size of this effect, but we believe it is not larger than 
-10%. 

Debye-Waller fits of the data of Table II were made to 
extrapolate back to the area at 0 = 0. For the Gaussian 
component it was found that 5 g (0) = 39.8 ± 7.0 and ( u h 2 ) 
= 0.37 ± 0.13 A 2 . (Note that (u b 2 ) characterizes the 
Gaussian contribution from the bound water and the 
agarose protons.) For the Lorentzian component, Lj, which 
is identified with mobile water, it was found that A m = 
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TABLE II 

AREA OF THE GAUSSIAN AND LORENTZIAN 
COMPONENTS IN ARBITRARY UNITS FOR A 20% 
agarose-h 2 o GEL 


Q 

Gaussian area 

Lorentzian area 

A-' 



0.7 

37 ± 5 

60 ± 19 

1.0 

27 ± 2 

64 ± 10 

1.2 

16 ± 1 

50 ± 7 

1.5 

19 ± 1 

37 ± 6 

1.9 

14 ± 1 

14 ± 5 


S m ( 0) = 81 ± 11 and (u m 2 ) = 0.39 ± 0.10 A 2 . In the 
deuterium oxide gel it was found that S( 0) = 24.4 ± 2.2 for 
the Gaussian component. Therefore, the area due to bound 
water, A b , is the difference between this area and the area 
of the total Gaussian in the water gel, and A b = 15.4 ± 7.3. 
The fraction of water that is bound is then f b = 0.160 ± 
0.066. Another way of expressing the amount of bound 
water is in terms of the number of grams of water bound 
per gram of solid (i.e., the bound hydration, c). This was 
found by multiplying the bound fraction with the hydration 
(4 g water/g agarose). The bound hydration is c = 0.64 ± 
0.26, which is in good agreement with the value c = 0.59 
quoted by Derbyshire and Duff (25). 

The mean-squared thermal displacement for the mobile 
water was the same as that obtained for pure water, within 
experimental uncertainty. Another useful comparison 
between the mobile water and pure water is the depen¬ 
dence of the width T of the Lorentzian on Q 2 , since this 
gives information about the diffusion parameters. This is 
shown in Fig. 3. The Lorentzian line width differs signifi¬ 
cantly from that of pure water for 0-values > 1 A 1 . The 
fact that the width for mobile water in agarose seems to be 
leveling off at a value significantly less than that of pure 
water indicates an increase of the residence time r (see Eq. 
3 ). 

Results for Brine Shrimp 

It is possible to hydrate the brine shrimp with deuterium 
oxide and maintain viability. However, it seems likely that 
substantial proton-deuteron exchange would occur in such 
a complex system, thus introducing uncertainties in the 
subtraction procedure. (Indeed, some exchange undoubt¬ 
edly occurs in the agarose gels [37].) For this reason we did 
not attempt to measure the deuterium-water difference 
spectra as we did for the agarose. Instead, the scattering 
from the brine shrimp was studied as a function of 
hydration. The data were then fitted with a single Gaus¬ 
sian representing the elastic scattering of both water and 
solid protons, plus a single Gaussian-Lorentzian convolu¬ 
tion representing the quasi-elastic scattering of both water 
and solid protons. From these data the contribution of the 
solids in the brine shrimp were estimated by extrapolating 
to zero hydration and then subtracting from the total 



Figure 4 Debye-Waller plot for the Lorentzian and the Gaussian 
components for the 1.2 g/g Artemia cysts. 

scattering. We identified the Gaussian component of the 
difference spectrum with bound water and the Lorentzian 
component with mobile water. 

The areas of the Gaussian and Lorentzian components 
as a function of 0 for brine shrimp for a hydration of 1.2 
g/g are plotted in Fig. 4. The other hydrations gave similar 
results. For the Lorentzian component, the area appeared 
to be relatively independent of 0. To estimate the area of 
this component, the area for the lowest 0 values were 
averaged together. The decision of how many 0 values to 
use was made on the basis of visual inspection of the data 
for each hydration. The number of points averaged on the 
log area vs. Q 2 plots were 5, 4, 3, and 3 for the hydrations 
0.10, 0.31, 0.78, and 1.2 g/g, respectively. 

For the Gaussian component, the areas seemed to follow 
a Debye-Waller dependence on 0 with an increase in slope 
as 0 —► 0. If the low 0 point was omitted, the areas could 
be fitted with a single Gaussian in 0. The percent of the 
QNS spectrum that was estimated to be in the Gaussian 
and Lorentzian components is shown in Table III. 

The percent of the QNS spectrum in the Gaussian is 
plotted against the concentration of cyst water in Fig. 5. A 
straight line fits these data well, and indicates that 86 ± 1% 
of the QNS spectrum would be Gaussian if there were no 
water present in the sample. This allows us to estimate the 

TABLE III 

DISTRIBUTION OF THE QNS SCATTERING 
BETWEEN THE GAUSSIAN AND THE 
LORENTZIAN AT 0 = 0 FOR THE ARTEMIA CYST 


Hydration 

% water 
by weight 

Gaussian 

Lorentzian 



% 

% 

0.10 

9.1 

78 ± 2 

21.6 ± 0.7 

0.31 

23.7 

65 ± 5 

35.3 ± 0.8 

0.78 

43.8 

48 ± 4 

52 ±4 

1.20 

54.6 

38 ± 3 

62 ± 5 


If the low Q point is included in the analysis, the area of the Gaussian is 
increased by -20%. 
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c 

0 



Figure 5 Percent area contained in the Gaussian component for 
Artemia cysts at Q - 0 vs. cyst hydration in grams per gram. Extrapola¬ 
tion to zero hydration gives 86% Gaussian contribution to the line area 
due to the solids. 


distribution between the Gaussian and Lorentzian compo¬ 
nents for the brine shrimp, assuming that this distribution 
did not change as a function of hydration. 

An elemental analysis of Artemia cysts shows (38) that 
for every gram of dry brine shrimp, there are 0.0801 grams 
of hydrogen. Because each gram of water contains 0.1119 
grams of hydrogen, we can estimate the amount of scatter¬ 
ing at each hydration that is due to the water and to the dry 
solids in the cysts. These estimates are shown in Table IV. 
The data are presented in terms of the percentage that is 
due to each component. 

Table III gives the percent of brine-shrimp solid 
scattering that is Gaussian and Table IV the percent of the 
QNS scattering that is due to the brine shrimp solids. From 
these we could calculate the amount of scattering to expect 
in the Gaussian and Lorentzian components due to the 
brine-shrimp solids. We substracted these figures from the 
data in Table III to calculate the amount of the Gaussian 
scattering due to water and the amount of Lorentzian 
scattering due to water. When we identified the Gaussian 
component of water with bound water and the Lorentzian 
with mobile water, we could calculate the bound fraction, 
as was done for the agarose gels. In addition, from the 
hydrations and bound fractions the values of c was calcu¬ 
lated. These results are given in Table V. 

This technique of analyzing the data did not allow us to 
compare the Q dependence of the line width of the mobile 
water Lorentzian with bulk water, as was done for agarose 
in Fig. 3. This was because we had not separated the 
contributions of the solids and the water at every Q value. 
We had only separated the contributions at Q = 0. We can, 
however, examine the Q dependence of the width of the 
Lorentzian component found for the highest hydration. 
The previous analysis leds to an estimate that at Q = 0 the 
solids contributed only -8% of the area of the Lorentzian. 
(37% of the scattering was due to the solids [Table IV], of 


TABLE IV 

PERCENT OF SCATTERING DUE TO EACH 
COMPONENT OF THE BRINE SHRIMP 


Hydration 

Solid scattering 

Water scattering 

g/g 

% 

% 

0.10 

87.74 

12.26 

0.31 

69.78 

30.22 

0.78 

47.85 

52.15 

1.20 

37.36 

62.64 


which 14% was in the Lorentzian line [Fig. 5]. The 
Lorentzian contribution due to the solids was thus 5% of 
the total, or 8% of the Lorentzian contribution of 62% from 
the water plus solids of the hydrated brine shrimp.) If we 
assume that the Debye-Waller factors of the Lorentzian 
and Gaussian terms did not differ greatly so that this 
contribution did not change significantly as Q increased, 
then the Lorentzian will be dominated by the mobile water 
at all Q values. Fig. 3 shows the width T(g) of the 
Lorentzian for h = 1.20, without any correction for the 
contribution to the Lorentzian from the bound protons. 
The width of this Lorentzian deviates dramatically from 
that of pure water, since it is relatively independent of Q 
and approximately three times smaller than that of pure 
water at large Q . The data point at Q = 0.5 A -1 was taken 
to confirm this behavior. The number of counts collected at 
this value of Q was increased to three times that collected 
at other Q values to reduce the statistical uncertainties. We 
did not take data at this low Q value for agarose or water, 
since the line became much narrower than the spectrome¬ 
ter resolution. However, the line from the Artemia cysts 
did not narrow. The line width vs. Q evidently cannot be 
well-described by the jump-diffusion model that is appro¬ 
priate for pure water. 

A determination of the Debye-Waller factor for the 
water in Artemia cysts at various hydrations was not 
possible due to the previously mentioned difficulty of 
subtracting the contribution from the solids at Q ¥= 0. For 
h = 1.2 g/g, we can obtain a mean-squared displacement 
for the “Gaussian component” from Fig. 4, but the areas 

table v 

PERCENT OF TOTAL SCATTERING DUE TO WATER 

COMPONENTS AND CALCULATED AMOUNT OF 

BOUND WATER IN THE BRINE SHRIMP 


Hydration 

Gaussian 

Lorentzian 

A 

c 


% 

% 



0.10 

2.8 ± 2.5 

9.5 ± 0.7 

0.23 ±0.16 

0.023 ± 0.016 

0.31 

4.5 ± 5.4 

25.7 ± 0.8 

0.15 ± 0.15 

0.046 ± 0.047 

0.78 

7.2 ± 4.0 

45 ± 4 

0.14 ± 0.07 

0.11 ±0.05 

1.20 

5.5 ± 2.5 

57 ± 5 

0.09 ± 0.04 

0.11 ±0.04 


Although the statistical errors are large, there is some indication that the 
bound hydration c varies with hydration, being smaller at the lower 
hydrations. 
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for this component contained a substantial (-50%) contri¬ 
bution from the protons in the solid as well as the contribu¬ 
tion from the bound water. We could not determine a value 
for the mean-squared displacement of the mobile water, 
since the rotational motion must be considered, as 
described in the next section. A comparison of the mean- 
squared displacements for water, 20% agarose, and 1.2 g/g 
brine shrimp is given in Table VI. 

DISCUSSION 

Rotational Contribution to the Scattering 

Law 


1.9 A -1 . This broadening was significant at high 0, but not 
much outside the error bars of our measurements. 

The same conclusions hold for agarose. Previous results 
(25) have shown that the NMR measured diffusion coeffi¬ 
cient Z> nmr in a 20% agarose gel is reduced to -0.62 Z> 0 , 
where D 0 is the value for pure water. A portion of this 
reduction can be ascribed to obstructions. If it is assumed 
that the polymers in the gel are equivalent to an ordered 
array of cylinders, then the obstruction effect can be 
estimated by (19, 39, and Rorschach and Hazlewood, 
unpublished data) 


D 


nmr 


D(agarose) 
1 + 0.80 


U ~ /b)> 


The results presented in Analysis and Results for the 
diffusive properties of pure water and 20% agarose gel are 
based on a particular model for the protonic motion: the 
jump-diffusion model as expressed in Eqs. 1 and 3. The 
parameters [areas, T(0 2 ) and the derived quantities D and 
r] obtained by fitting the experimental data to this model 
by the least-squares method, therefore cannot be regarded 
as really fundamental, since they are model dependent. 
Nevertheless, we believe that the comparisons made 
between the systems have validity beyond this particular 
model, since we have consistently applied the same model 
in the same way. 

In the analysis of our data, we used only a single 
Lorentzian to fit the scattering law given by Eq. 5. This 
scattering law consists of two Lorentzian lines L t and L r+t 
of different widths. In the case of water, the influence of 
the second line, L r+t , which was broadened by the rota¬ 
tional motion, lead to a small increase in the width 
obtained by a least-squares fit of the data with a single 
Lorentzian. At low 0, where the translational line was 
narrow, the relative amplitude of the rotationally broad¬ 
ened line was small, while at high 0, the translational line 
was very broad, so that the widths of the two lines were 
more nearly the same. Calculations based on this model 
with a value of D r = 8 x 10 10 s _1 for pure water (see 
Discussion section on Comparison with NMR) showed 
that the net result was that the measured line was broad¬ 
ened by -4% at 0 = 0.5 A -1 , increasing to -15% at 0 = 

TABLE VI 

MEAN-SQUARED DISPLACEMENTS FOR THE 
GAUSSIAN (ELASTIC) AND LORENTZIAN 
(QUASI-ELASTIC) CONTRIBUTIONS FROM THE 
LISTED SCATTERERS 



<w b 2 > 

elastic 

<U m 2 > 

quasi-elastic 


A 2 

A 2 

Pure Water 

— 

0.33 ± 0.02 

H 2 0 in 4 g/g agarose 

0.37 ± 0.3 

0.39 ± 0.10 

H 2 0 in 1.2 g/g brine shrimp 

0.23 ± 0.03* 

— 


*Contains contributions from solids. 


where Z>(agarose) is the “local” diffusion coefficient in 
agarose, <f> is the volume fraction of the cylinders (solids), 
and / B = bound fraction = 0.16. For </> = 0.2, we find 
Z>(agarose) = 0.8 D 0 . A calculation of the influence of the 
rotational motion on the agarose line width can now be 
made, as we did above for pure water. If we assume a 
rotational diffusion coefficient equal to that of pure water, 
we find a broadening of the measured line by -4-20%, 
depending on the value of 0. This calculation probably 
overestimates the broadening, since it is likely that the 
rotational diffusion coefficient is considerably less in the 
agarose gel, which would greatly reduce the broadening 
caused by the rotational contribution. This will be consid¬ 
ered further in the Discussion section on Comparison with 
NMR. 

The Artemia results, in contrast, can only be understood 
by taking the rotational motion into account, since the 
translational diffusion coefficient was greatly reduced 
from that of pure water. Previous NMR measurements in 
Artemia (39) show that the diffusion coefficient was 
reduced by a factor of 7 for a hydration of 1.2 g/g. Some of 
this reduction may be due to the presence of diffusion 
barriers (e.g., membranes and macromolecular struc¬ 
tures). Indeed, the pulsed-field-gradient measurements of 
Tanner (40,41) show that approximately one-half of the 
reduction was due to barriers, leaving an estimated reduc¬ 
tion of 3.5 for the “local” diffusion coefficient. Such a 
greatly reduced value for D will have a strong influence on 
the observed neutron scattering. In fact, the translational 
line will now be much narrower than the spectrometer 
resolution, and it will be included with what we have 
termed bound water in our analysis. The observed Lorent¬ 
zian line will be the rotationally broadened term (the 
second term of Eq. 5), and the observed width will be r r + 
T t . The diffusion parameters obtained from this analysis, 
based on Eq. 5, are presented in Table VII. The curves 
drawn through the data points in Fig. 3 are fits with these 
parameters. 

The 0-dependence of the area of the Lorentzian line 
(see Fig. 4) is not inconsistent with this analysis. According 
to Eq. 5, the area should be proportional to j 2 (Qa) exp 
(— 0 2 ( t/ m 2 )). This factor approaches zero as 0 —► 0, but it 
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TABLE VII 

DIFFUSION PARAMETERS USED TO FIT 
SCATTERING LAW TO QNS DATA* 


System 

D/D 0 

V T o 

D, 

Pure H 2 0 

1 

1 

io’s-' 

20% agarose gel (4 g/g) 

0.80 

1.4 

— 

Artemia cyst (1.2 g/g) 

0.29 

3.9 

6.0 


*See Fig. 3. 


has a very broad maximum at Q « 1 A -1 for (u m 2 )/a 2 = 

I /2 and a =* 1 A, and is relatively flat in the entire Q range 
of our experiments. 

Multiple Scattering Corrections 

A possible source of error in relating the line widths to the 
diffusive parameters through Eq. 3 is multiple scattering 
effects. Multiple scattering will alter the Lorentzian line- 
shape predicted by Eqs. 1 and 3, and least-squares fits to 
the line by a single Lorentzian will in general lead to a 
broader line than would be obtained from single scattering. 
The line area from which the Debye-Waller factor and 
hydration fractions are determined will also be affected. 
Multiple scattering effects have usually been ignored in 
quasi-elastic scattering experiments. This has been justi¬ 
fied on the basis that measurements with different sample 
thicknesses give similar results, or that the effects are 
negligible in “sufficiently” thin samples (e.g., with trans¬ 
missions -90%), or that inelastic background subtraction 
compensates for it. Such arguments are dangerous, since 
the influence of multiple scattering depends on sample 
shape as well as on thickness, and the influence of sample 
shape on line areas may be appreciable. 

The best way to determine the influence of multiple 
scattering is to calculate the scattering with a Monte-Carlo 
method. J. R. D. Copley has developed a program called 
MSCAT (42,43) that calculates the observed scattering 
based on a theoretical scattering law and a sample geome¬ 
try supplied by the user. We have used this program to 
determine the scattering due to the quasi-elastic portion of 
the spectrum for an experimental geometry similar to that 
used in this study. The scattering law described by Eqs. 1 
and 3 was used with D = 2.4 x 10“ 5 cm 2 /s and r — 10 -12 s. 
The “sample” used in the program was a 1-mm thick 
slab-shaped sample of pure water. Simulated data were 
obtained for this “sample” from the MSCAT calculation. 
These data were treated as experimental data, and our 
fitting routines were used to determine the effective areas 
and widths. 

Fig. 6 illustrates the corrections to the width and 
integrated intensity which these calculations predict. The 
errors are on the order of 10% in the width and do not vary 
much with Q . The amplitude corrections are of the same 
order, except near Q = 1.5 A -1 , where they are considera¬ 
bly larger. There is no strong experimental evidence for the 
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Figure 6 Multiple scattering correction for a 1-mm slab of H z O as 
determined by a Monte Carlo calculation: ( a ) Ratio of line area with 
multiple scattering to that without multiple scattering vs. Q. (b) Ratio of 
line width with multiple scattering to that without multiple scattering vs. 

Q 

large decrease in area at Q = 1.5 A -1 (see Fig. 2), and we 
have made no corrections to any of our data for multiple 
scattering. Because all samples had nearly the same trans¬ 
mission, the correction would be roughly the same in all 
cases, so that, although the absolute values of the parame¬ 
ters derived from the widths and areas may be in error by 
-10%, the relative values should be unaffected. 

Comparison with NMR and Dielectric 
Relaxation Results 

The present neutron scattering results combined with 
previous NMR data on agarose gels and Artemia cysts 
give new insight into the disposition and dynamics of the 
water in these systems. 

Agarose Gels. The neutron scattering can be 
well described by a jump-diffusion model in which D — 0.8 
D OJ t = 1.4 r 0 ( D 0 and r 0 are the values for pure water.) The 
line plotted in Fig. 3 was obtained from Eq. 3 with these 
parameters. This value for D is consistent with the NMR 
results on diffusion (25) and indicates that there is a 
“local” reduction in the values of D that cannot be 
associated with obstructions or barriers. This might not be 
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evident from the initial slope of T(0 2 ) in Fig. 3, but that 
conclusion is justified by the increased value of the resi¬ 
dence time r that is related to the diffusion coefficient by 
6Dt = (/ 2 ), where (l 2 ) is the mean-squared jump 
distance. The influence of the gel is thus to stablize the 
water structure without much change in average water- 
water distances but with an increased “binding” time 
between jumps. 

The mean-square atomic displacement did not differ 
greatly from that of pure water (Table VI), which indi¬ 
cates that, although the residence time was increased due 
to association with the polymer chains, the local oscillatory 
motion, and therefore the binding forces experienced by 
the protons, were not much changed. 

The NMR relaxation time T x has been measured by 
Woessner and Snowden (37) for agar gels and by Derby¬ 
shire and Duff (25) for agarose gels. Woessner and 
Snowden show that (1/T,) agar = (l/r o )° + 8.2 Cat 25°C, 
where C is agar concentration in grams per cubic centime¬ 
ter of H 2 0, and the zero superscript indicates the value for 
pure water. This gives a T x of 0.5 s for a 20% gel, in good 
agreement with an extrapolation of the reesults of Derby¬ 
shire and Duff, which gives T x ^ 0.4 s for a 20% agarose 
gel, a reduction of a factor of -6 from the value for pure 
water. 

The T x value for the water in agarose gel contains 
contributions from both rotational and translational 

motion (23). (1 / -Fl) a garose = 0/^l)rotational 4* (1 / ) translational* 

The rotational relaxation rate was determined by the 
rotational diffusion coefficient D x and the translational 
relaxation rate by D. In the shortcorrelation time limit 
(23), (1 / 7",)i a 1 /D x . Thus, the 20% reduction in D will not 
produce the observed reduction in T x . The value of D r for 
the gel is unknown, but the analysis of the Artemia data in 
the next paragraph suggests that D r may be considerably 
smaller in the gel than in pure water. If D r were reduced by 
a factor of -8, then the value of T x = 0.5 s could be 
explained as a bulk mechanism without invoking rapid 
exchange with a solidlike hydration layer. Such a reduction 
in D r would also account for the good fit of the neutron 
scattering data to the jump diffusion model, since the 
rotational motion would then broaden the line by only a 
few percent. In this picture, T x would be due to bulk 
mechanisms, while T 2 would be determined by exchange 
with a solidlike phase whose correlation time is too long to 
affect T x . 

Artemia Cysts. The neutron scattering in this 
case can be well described by a combination of transla¬ 
tional and rotational diffusion for which D = 0.29 Z> 0 , r = 
3.9 r 0 , and D r = 6 x 10 9 s" 1 . The line plotted in Fig. 3 is 
F((? 2 ) = T r + r t , where T r = 4 D T and T t is given by Eq. 3, 
with these parameter values. This value of D is consistent 
with the NMR diffusion measurements (39), and shows 
that the large reduction in D is not entirely due to barriers 


or compartments, but is a true local effect due to the 
stabilizing of the water by the nonaqueous components 
present in these cells. The picture is very much the same as 
for agarose, except more pronounced due to the relatively 
higher concentration of solids in the Artemia system and a 
consequent larger reduction in D. The amount of intercel¬ 
lular surface area in living cells is indeed very large (3), so 
that the result is not surprising. The mean-square atomic 
displacement is considerably less than that of pure water, 
but this may be due to the large (-50%) contribution from 
the bound protons (see Table VI). 

The value of D r for Artemia can be compared with 
results from dielectric relaxation studies by Clegg et al. 
(44). The comparison is not entirely straightforward, since 
dielectric dispersion measurements do not measure D r 
directly. Furthermore, the problem of exchange among the 
various fractions of water also needs careful consideration 
in interpretations of dielectric results. 

Previous measurements on pure water (45) have shown 
that the Debye relaxation time r D is 0.93 x 10“ 11 s at 20°C. 
The “microscopic” relaxation time is smaller, due to 
local field effects, and is (reference 45, p. 299) r M =* r D /1.5. 
The rotational diffusion coefficient is related to r M by (23) 
D r = 1/2 r M , and thus D r = 8.1 x 10 10 s -1 for pure water. 

The value Z) r = 6 x 10 9 s" 1 for Artemia measured by 
neutron scattering is lower than that of pure water by a 
factor of -13. A similar increase in r D has been observed in 
systems of relatively low water content, such as hemoglo¬ 
bin (46, 47) and lysozyme (48). These changes are much 
larger than those obtained by Clegg et al. (44), who have 
measured the dielectric dispersion for water in Artemia 
cysts for frequencies up to 70 GHz. They find a relaxation 
time only slightly longer (-2 times) than that of pure 
water, which suggests D r is reduced by a similar factor. 
Previous work on biological systems of higher water con¬ 
tent (49) also reported values of r D close to that of pure 
water. Nevertheless, this is not convincing evidence that 
the water in the cell is predominantly “free” water. Should 
rapid exchange occur, the dielectric properties would be 
dominated by the presence of only a small fraction of free 
water, in the same way that the value of T 2 in biological 
systems is influenced by a small rapidly relaxing phase of 
water. Such considerations have not been taken into 
account in the analysis of microwave dielectric studies of 
water in cells and tissues (44,49), and they may be 
important. 

The above values for D and D r obtained from the 
neutron spectra are also consistent with the NMR relaxa¬ 
tion time T x . The value of T x for cysts at 1.2 g/g water 
content is -250 ms (50), a factor of -12 less than that of 
pure water. For pure water, T x is determined by the 
rotational and translational motion (23) 

fF - (fF ♦ (fF • 

1 1/ ' 1 1 / rotational V 1 1 / translational 
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where (51) 


rotational V -* 1 

The neutron results show that the rotational contribution 
to (1/7,) for Artemia cysts was -13 times that for pure 
water, since the value of D x was reduced by this factor. The 
translational contribution was increased by a factor of -4. 
This gave an expected reduction factor in T x for Artemia of 
-10, which agrees rather well with the observed values of 
7, (50). 


/ translational 



CONCLUSIONS 

We thus conclude that the neutron, NMR, and dielectric 
relaxation measurements on Artemia cysts are all consis¬ 
tent with a picture in which the majority of the cell water 
has strongly reduced translational and rotational diffusion 
coefficients that are not due to obstructions, compart¬ 
ments, or exchange with minor phases, but are instead an 
intrinsic feature of the intracellular water. This result 
supports hypotheses, such as those of Berendson (52) and 
Ling (5-10) that the motional freedom of water molecules 
is reduced by association with cellular macromolecules. 
The consequences of this conclusion could be far reaching, 
since much of current thought about cell structure and 
function considers the water to exhibit properties that are 
the same as those of bulk aqueous solutions. It seems very 
likely that the role of water in chemical reactions, protein 
conformations, ion exclusion, and other cellular functions 
will prove to be a more dynamic one than has heretofore 
been assumed (5-10, 53). 
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